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Abstract 

We suggest a method to prepare any chosen superposition ao|0) + ai|l) of 
the vacuum and one-photon states. The method is based on a conditional 
double-interferometer fed by an one-photon state and a coherent state. The 
scheme involves only linear optical elements and avalanche photodetectors, 
and therefore it should be realizable with current technology. A realistic 
description of the triggering photodetectors is employed, i.e. we assume that 
they can only check, with a certain efficiency, whether or not any photon 
is present. We discuss two working regimes, and show that output states 
with fidelity arbitrarily close to unit may be obtained, with non vanishing 
conditional probability, also for low quantum efficiency at the photodetectors. 

I. INTRODUCTION 

The last two decades have witnessed a substantial development in quantum engineering 
and measurement of light. Several kinds of nonclassical states of light can now be generated, 
and their quantum properties can be fully characterized by accessible measurement schemes 
W0|. Besides fundamental interest, nonclassical states also find applications, as for exam- 
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pie, the use of number states in quantum communication channels, and of squeezed light 
in high-precision gyroscopes and interferometers. More recently, the quantum engineering 
of light received new attention, which is mainly motivated by the potential improvement 
offered by quantum mechanics to the manipulation and the transmission of information || . 
Indeed, phenomena like teleportation ||, and quantum dense coding || found their first 
implementation in the quantum optical domain. 

Photons do not interact, and this feature is very useful for the transmission of infor- 
mation without signal degradation. Indeed, the typical figures for losses in optical fibers 
are below 0.3 dB/Km. On the other hand, the same characteristic poses limitations to the 
manipulation of the quantum information encoded into a quantum state of light. Photon- 
photon interactions needed for computation, in fact, take place only in active optical media, 
characterized by nonlinear susceptibility. Usually, such nonlinearities are small, or masked 
by the concurrent absorption processes. Only recently, new methods based on dark atomic 
resonance and electromagnetically induced transparency || have been suggested to strongly 
enhance nonlinearity while suppressing absorption. The possibility of such giant nonlinear- 
ities renewed the interest for optical quantum technology, as it opens new perspectives to 
build single-photon quantum logic gates. 

In this paper, we devote our attention to the preparation of any chosen superposition of 
the vacuum and the one-photons states a |0) + ai|l). This is the simplest state of light that 
carries a complete phase information, and, in turn, it represents the simplest example of an 
optical qubit. Remarkably, this is a low-energy-expense encoding of quantum information, as 
it requires, in average, less than one photon for each qubit. In particular, for the conventional 
computational basis \i/j±) = ^75 [|0) ± |1)], i.e. for balanced superposition, half a photon for 
each qubit is required. 

Different methods have been discussed with the purpose of engineering superpositions of 
radiation states. Mostly, these are in the context of cavity QED, since the interaction with 
atoms passing through the cavity allows to select specific components of an initial signal [|7| . 
More recently, an all-optical device, based a ring cavity coupled to the signal through a Kerr- 
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medium, has been suggested to realize Fock filtering, and thus preparation of superpositions 
||. In addition, a conditional scheme based on beam splitters and photocounters has been 
suggested to implement the optical state-truncation of a coherent state, which, in turn, 
is used to prepare superpositions. As we will see, this last setup corresponds to a particular 
case of the present proposal. 

The present scheme involves only linear optical elements and avalanche photodetectors, 
and therefore it should be realizable with current technology. In essence, it consists of two 
Mach-Zehnder interferometers arranged such that one of the outputs from the first one is 
then used as input for the second one. The first MZ is fed by an one-photon state, whereas 
the second is fed by a weak coherent state. The output states of the second interferometers 
are then measured, and the conditional output state from the first MZ turns out to be a 
superposition of the vacuum and one-photon states. The amplitudes for the two compo- 
nents can be tuned by varying the internal phase-shifts of the two interferometers and the 
amplitude of the coherent input. 

The paper is structured as follows. In the next Section the scheme is presented and its 
dynamics is evaluated. Then, the ideal conditional output state is calculated for the case 
of a perfect photodetection process. In Section flTJ we take into account the imperfections 
of realistic photodetectors, and study their effects on the preparation of the superposition. 
In the literature two models of photodetectors (PDs) have been employed. In the first, 
which we use throughout the paper, it is assumed that PDs are only able to check, with 
a certain efficiency, whether or not photons are present. This is a reliable description of 
customary avalanche photodetectors, and we refer to this as YES/NO photodetection. In 
the second model, PDs are still affected by non unit quantum efficiency, but now they are 
able to discriminate among the number of incoming photons, thus acting as photocounters. 
This description does not yet correspond to available PDs, as no evidence of detectors 
capable to discriminate between the presence of, say, n and n + 1 photons for a generic n, 
have been reported. The photocounter model has been sometimes used in the literature 
on conditional measurements, remarkably in Ref. ||, and this led to the conclusion that 
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schemes are reasonably insensitive to the detectors' inefficiency. In general, this is no longer 
true when the realistic features of avalanche photodetectors are taken into account. On the 
other hand, we will show that the present scheme offers a working regime in which the use of 
YES/NO detectors is enough to assure the reliable preparation of any chosen superposition. 
Finally, Section [TV] closes the paper with some concluding remarks. 

II. CONDITIONAL DOUBLE-INTERFEROMETRY 

The scheme we have in mind is the conditional double-interferometer depicted in Fig. 
|TJ. It consists of two Mach-Zehnder interferometers in cascade, in a way that makes one of 
the output signals from the first one to constitute one of the input signals for the second 
one. The three field-mode involved in the setup are denoted by a, b, and c, whereas the 
BS's are symmetric beam splitters. We also assume that equal and opposite phase-shifts, 
denoted by 9\ and 62 respectively, are imposed in the arms of each interferometer (see Fig. 
[l]). We assume that two interferometers are built with identical balanced beam splitters, 
and this means that they are fully characterized by the value of the internal phase-shift 
between their arms. The first interferometer is fed by an one-photon state in the mode b, 
whereas the other port, corresponding to mode a, is left unexcited. After this first stage the 
photon has a nonzero amplitude of being in both the output paths, whereas the values of 
such amplitudes are determined by the internal phase-shift 9\. The output mode b is then 
mixed with mode c, prepared in a weak coherent state I7), and they are both detected at 
the output of the second interferometer. Depending on the result of the measurements we 
have different conditional output states for the output mode a of the first interferometer. 
In particular, we will see that any chosen superposition ao|0) + ai|l) of the vacuum and 
one-photon states can be prepared with nonzero conditional probability. 

The evolution operator of each interferometer can be written as 

Vmz(9) = Vbs exp {iO(a!a - b ] b)} V ] BS , (1) 

where 
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Vksr = exp{z^(a t 6 + & t a)} , (2) 

denotes the evolution operator of a symmetric beam splitter. Eq. ([!]) can be written as 

Vmz{4>) = ex P j^^&j exp {— i(f>(atb + tfa)} exp j— *7^&j > (3) 

which shows that a Mach-Zehnder interferometer is equivalent to a single beam splitter BS^ 
of transmissivity r = sin 2 0, where = 9/2, preceded and followed by rotations of n/2 



performed on one of the two modes flO ]. After straightforward algebra, one may write the 



evolution operator of the whole device as 

U(4>i,4>2) = exp j«— b^bj exp {i(j) 2 (tfc + c*&)} exp {i<pi(a)b + &*a)} exp |— i— M&j . (4) 
The overall input state can be written as 

|* IN > = |0) a |l> 6 |7> 6 , (5) 
whereas using Eq. (f|) we obtain the expression for the overall output 

I^out) = C/(0i,0 2 )|*i N ) = cos 0i |l) a |7cos0 2 ) fe |7sin0 2 ) c 

+ sin 0i sin 02 |0) o 6 t |7cos0 2 )b|7sin0 2 ) c 

-sin 0i cos 2 |O) a |7cos0 2 ) 6 c t |7sin0 2 ) c . (6) 

We are now ready to analyze the effect on the mode a of a measurement performed on the 
output modes b and c. For the moment let us assume that the generic measurement of a 
quantity X on the b mode, and Y on the c mode is performed; we also assume that the two 
quantities are independent on each other. The measurement is thus described by a factorized 
probability operator measure (POM) fl = Il x ® fl y , where x and y denote the possible 
outcomes for the two quantities. The POMs are positive (hence selfadjoint) operators, if 
X and Y denote the spaces of the possible outcomes the normalization conditions can be 
written as (I denotes the identity operator) 

dxfl x = I I dy fly = I 



The probability of the event (x, y) is given by the global trace over the three modes 



p tv 

1 xy - LL abc 



I^outX^outI n* <8» ft y ] , (7) 



whereas the corresponding conditional output state for the mode a is given by the partial 
trace 

Qxy = -p—Ttbc |*out)(*out| ti x <g> fl y . (8) 

xy L - 1 

Actually, different measurements on b and c lead to different conditional output states, and 
since some events are more likely to occur than other, so are the corresponding conditional 
output states. In the following we consider the measurement of the photon number, and 
in particular we focus our attention to the case of a single photon recorded in one of the 
output modes, and no photons recorded in the other one. In order to establish notation we 
assume that the photon is recorded in the output mode b. However, the results are valid 
also for the case of of a photon registered in the mode c, up to a tt/2 shift in the internal 
phase-shifts of the second interferometer. 

The ideal measurement of the photon number on the two output modes is described by 
the POM 

n„fe = \n) bb (n\ <g) \k) cc (k\ n, k = 0, 1, ... . (9) 
Therefore, the detection probability for the case n = 1, k = is equal to 

Pio = I (^out|1>6|0) c | 2 = e^l 2 [sin 2 fa sin 2 2 + | 7 | 2 cos 2 fa cos 2 fa] , (10) 



and the conditional output state for the mode a is given by 

OUT I |1)m,(1|®|0> cc (0| 



f?io — -5— Tr 6c 
-n.0 



Pio 

where 



= 4 [coo|0)(0| + Cll |l)(l| + c 01 |0)(l| + c 10 |l)(0|] , (11) 



c oo = e ' 7 ' sin 2 0! sin 2 02 C\\ = e ' 7 ' |7| 2 cos 2 0i cos 2 2 (12) 
c i = e~' 7 ' 2 7sin0! sin0 2 cos0i cos02 c w = c* Ql (13) 

the star denoting complex conjugation. By looking at Eqs. ([□]) it is easy to recognize that 
£10 is actually a pure state £10 = iV'ioXV'iol where 



|^ ^ _ sin0isin0 2 |O) + 7cos0i cos0 2 |l) ^ 
a/ sin 2 0i sin 2 2 + M 2 cos 2 0i cos 2 02 
By varying the interferometric shifts, and the amplitude 7 of the coherent input c, we may 

achieve any chosen superposition of the vacuum and the one-photon states. In particular, 

the internal phase-shifts 0i and 02, and the modulus I7I, govern the weights of the two 

component, whereas the relative phase equals the argument of the complex amplitude 7. 

We notice that the truncation scheme of Ref. || is equivalent to a particular case of our 

setup, corresponding to the balanced choice 0i = 02 = it /A. This is due, as mentioned above, 

to the fact that a Mach-Zehnder interferometer is substantially equivalent to a single beam 

splitter of transmissivity r = sin 2 0. However, besides the fact that the present scheme offers 

additional degrees of freedom, the use of an interferometric setup has the specific advantage 

of a larger stability. 

Remarkably, the superposition of Eq. (|H| ) may be obtained by different values of 0i, 
02 and 7, and this degree of freedom can be used to maximize the corresponding detection 
probability P10, i.e. the probability of the event which leaves the state of the mode a into the 
desired superposition. Let us consider, for example, the preparation of a generic balanced 
superposition 

m - • (15) 

In order to reduce Eq. ( JUty to Eq. flT5] ) we need if = arg7 and I7I = tan0itan0 2 . The 
detection probability is then given by 

P^ = 2 sin 2 0! sin 2 2 exp [ — tan 2 tan 2 ] . (16) 

From Eq. flTEp we have that < P^q < 0.21, with the maximum value reached for 0i = 02 ~ 
0.715, corresponding to an optimum amplitude |7| Q pt — 0.755. However, the dependence 



of the detection probability on 0i and 02 is not dramatic, such that there exists a sizeable 
region in which the detection probability is above Pj* > 20% (see Fig. |2|). 

As mentioned above, the symmetric case of a photon detected in the output mode c and 
no photons in b leads to an equivalent result, up to the replacement 02 — » 02 + 7f/2. In 
formula we have 

Pox = e~' 7 ' 2 [sin 2 X cos 2 2 + |7| 2 cos 2 0i sin 2 2 ] , (17) 

and 



|^ _ sin0! cos0 2 |O) -7COS0! sin0 2 |l) 
a/ sin 2 0x cos 2 02 + I7I 2 cos 2 0i sin 2 2 
By comparing Eqs. ( |I4D and (|T8| ) we also note that the scalar product between the two 

conditional output state is given by 



(■^oil^io) oc sin0 2 cos0 2 ^sin 2 0i - |7| 2 cos 2 0ij , (19) 

which means that for I7I = tan0i, and independently on 02, the two states are orthogonal 
(this also happens for 2 = pvr/2, p G Z and any 7 and 0i, but this case just corresponds 
to have one state in the vacuum and the other in the one-photon state). In this regime, 
the scheme provides a reliable source (i.e. with P 10 = P01 = exp{— tan0 2 } sin 2 X which is 
larger than 20% in the region around 0i ~ 0.67) of a quantum-optical computational basis. 

Of course, it is now of interest to study whether or not the superpositions of Eqs. (|I~4D 
and ( |18D may be obtained in a realistic implementation of the setup. Since the scheme is 
based on conditional measurements, the main concern should be with the photodetection 
process. Therefore, in the next Section we are going to take into account the imperfections of 
available photodetectors, in order to check the robustness of the preparation scheme against 
the detectors inefficiency 

III. EFFECTS OF REALISTIC PHOTODETECTION 



Light is revealed by exploiting the interaction with atoms or molecules. Each photon 
ionizes a single atom, and the resulting charge is then amplified to produce a measurable 
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pulse. In practice, however, available photodetectors are hardly performing the ideal mea- 
surement of the photon number. Their performances, in fact, are limited by two main kinds 
of imperfections. On one hand, photodetectors are usually characterized by a quantum ef- 
ficiency lower than unit, which means that only a fraction of the incoming photons lead to 
an electric pulse, and ultimately to a "count". Some photons are either reflected from the 
surface of the detector, or are absorbed without being transformed into electric pulses. On 
the other hand, customary photodetectors involve an avalanche process to transform a single 
ionization event into a recordable pulse. This implies that it is very difficult to discriminate 
between the presence of a single photon or more than one. 

The outcomes from such a detector may be either YES, which means a "click", cor- 
responding to any number of photons, or NO, which means that no photons have been 
recorded. This kind measurement is described by a two-value POM 



yi n = J2^-v) p \p)(p\ n Y = i-n r 

p=0 



(20) 



where rj is the quantum efficiency, and I denotes the identity operator. Indeed, for high 
quantum efficiency (close to unit) Hn approaches the projection operator onto the vacuum 
state, and fly onto the orthogonal subspace. 

The event of observing a click at the PD surveying the output mode b (i.e. Db, see Fig. 
lip, and no photons at D c , is characterized by the probability 



P YN [?7,7,0i,02] = Tr abc |^ ut)(^out| IIy ®n 



- V \y\' sin 2 4>2\ l _ e -nh\ 2 cos 2 <fe + ^ ^2 ^ 



+ COS <p2(j]\l\ sin 02-1 

The corresponding conditional output state is 

1 



-»?|7| 2 cos 2 02 



R 



-Tr 



be 



YN 
1 



I^outK^outI n Y ®n P 



Py 



doo |0> <0| + d u |1)(1| + d i |0)(1| + d* m |1)(0| 



(21) 



(22) 



where the coefficients are given by (see appendix A) 



du = e 
d o = e 



-??|7| 2 sin 2 02 



-?7|7| 2 sm 2 02 01 -„2 



COS (pi 

sin 2 0i 



1 - e 



-r;|7| 2 cos 2 02 



1 - (1 - r])e-^ 2 cos2 * 2 + r] cos 2 2 (??|7| 2 sin 2 2 - 1) 



c?oi = e ^l 7 ' sm ^ 777 sin (pi sin 2 cos (pi cos 2 



(23) 



In general, the conditional output state g YN is no longer a pure state. However, as we will see, 
there are regimes in which g YN approaches the desired superposition. In order to compare 
Qyn with the ideal conditional output \ipio) we consider the fidelity F = (V'ioI^ynIV'io)- From 
Eqs. (D and (||423]) we have 



F[V,1, 



1 



P Y n sin 2 (pi sin 2 (p 2 + I7I 2 cos 2 (pi cos 2 (p 2 
+2?]|7| 2 sin 2 0! sin 2 (p 2 cos 2 (pi cos 2 (p 2 + sin 4 0! sin 2 (p 2 

+?7cos 2 2 (7y|7| 2 sin 2 2 - 1 



2 .... 1 :.,2 , / 1 „-»?|7| 2 cos0| 



cos ' o :[ sm" (p 2 \l — e 

1 - (1 -^e^H 2 ^ 2 ^ 2 
(24) 

Our goal is now to find regimes in which the fidelity of the conditional output state is 
close to unit and, at the same time, the corresponding detection probability P YN does not 
vanish. In particular, we are interested in the preparation of those superpositions where the 
amplitudes of the two components are of the same order, thus assuring that the state is far 
from being just the vacuum or the one-photon state. This requirement roughly corresponds 
to the condition 



sm (pi sm (p2 — cos (pi cos (p 2 



(25) 



whereas a fine tuning of the amplitude, as well as the phase of the superposition, may be 
obtained by varying the complex amplitude 7 of the coherent input. The condition in Eq. 
(pop is satisfied by two different working regimes of the setup, i. e. by two different pairs 
of values of the internal phase-shifts. These are the balanced choice (pi = (p 2 = 7r/4 and the 
unbalanced one (pi ~ (p 2 — tt/2 — (pi respectively. In both cases, the general expression 
for the fidelity in Eq. fl2~4] ) may be considerably simplified, and the corresponding working 
regime discussed with some details. 
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A. The case 4>i = 4>2 = tt/4 



In the case 4>i = <p2 = tt/4, the detection probability rewrites as 

P YN [ V , 7, vr/4,vr/4] = e"^ 2 jl - + -^[e^ 2 + ^| 7 | 2 - 2)] | , (26) 

and the fidelity 

A-2t]+ | 7 | 2 (2 + 7?) 2 -Ae~^ 2 (l - rj + | 7 | 2 ) 
Ffo, 7, tt/4, tt/4] = (1 - | 7)2) [g - _ 2) _ 4e -,| 7 p/2( 2 _ ^ • ( 27 ) 

The detection probability shows relatively large values (P YN > 50%) in a sizeable region 
(see Fig. |3|) of the 7] — | 7 | space, including also situations with low quantum efficiency. 
Unfortunately, the fidelity of Eq. ( p7|) is a rapidly decreasing function of the coherent 
amplitude and, in the relevant region < I7I 2 < 4, it is bounded by 

F c 2 + M 2 (9-4e-HV 2 ) 

^ MAX< (l + | 7 |2 )( | 7 |2 + 6 _ 4e -| 7 |V2 ) • W 

This working regime is thus effective only for | 7 | <C 1, corresponding to the preparation of 
superposition where the vacuum component is preponderant. Indeed, for balanced superpo- 
sitions (| 7 | ~ 1) we have the bound F < 93%. As we will see in the following, this limitation 
can be overcome by the unbalanced tuning of the internal phase-shifts. 

Our analysis of the balanced scheme led to conclusions that are in contrast with those 
of Ref. 0, where, as mentioned in Section |I|, a formally equivalent scheme has been used. 
The reason for this disagreement stays in the different models employed to describe the 
photodetection process (see Appendix |B|). Actually, as far as we know, the YES/NO model 
used here is more realistic that the photocounter model used there, as, in fact, no evidence 
of detectors capable to discriminate between the number of incoming photons have been 
reported. We should conclude that the authors' hope of "reasonable insensitivity" to the 
detectors inefficiency |J is not yet realized with current technology. 

B. The case (f>i ~ , <p2 — k/2 — <j)\ 

In the case <j>\ ~ 0, <p2 — tt/2 — <fii, the detection probability is given by 
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vh\ 2 



P YN [?7, 7, 0i ~ 0, 02 ^ vr/2 - 0i] ~ e 



^7 1 | 2 cos 2 02 + i] sin 2 0i 





' -vh\ 2 


7 2 ) exp 





(29) 



and the fidelity 

F[ri, 7, 0i ^ 0, 02 ~ vr/2 - X ] ~ 1 - ^0 2 (2 - r,) . (30) 

Remarkably, the fidelity is now independent on the amplitude of the coherent input, and 
can be made arbitrarily close to unit by choosing a smaller value for 0i. The price to pay 
for this result is a lower value of the detection probability, namely a lower efficiency of the 
preparation scheme. However, the resulting probability is still large enough to make the 
scheme an effective source of superposition states. As an example, let us consider 0i such 
that sin 2 0i = cos 2 02 — 0.01. In this case, we obtain a very high value of the fidelity 
F > 99%, and yet a detection probability given by P YN ~ 1% (almost independently on the 
quantum efficiency). More generally, we can substitute Eq. ( f29|) in Eq. (|3~0|) to write 

Eqs. (p9"D , ( |30"D and ( plf ) assure that a reliable generation of the desired superposition is 
achievable (with non vanishing probability) also for low quantum efficiency at photodetec- 
tors. 



C. Balanced superpositions 

We end the Section by illustrating the performance of the setup in preparing the special 
class of exactly balanced superpositions of Eq. (|1^). The requirement for equal amplitudes 
reads |7| = tan 0i tan 02. By substitution in Eqs. (pl|) and ( p2|) we obtain the detection 
probability and the corresponding conditional output state. We do not show here the result- 
ing expressions, which are rather cumbersome. Instead, in Fig. [|, we report the behavior of 
both, the fidelity and the detection probability, as a function of the internal phase-shifts for 
two values of the quantum efficiency at the photodetectors. As it is apparent from the plots, 
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there always exists a region, in which the fidelity is very close to unit, and yet the detection 
probability is larger than 10%. Therefore, for balanced superpositions, the performances 
with realistic detectors do not substantially differ from that obtained in the ideal working 
regime discussed in Section [TI|. 

IV. CONCLUSIONS 

We have analyzed a linear, conditional, interferometric setup to prepare any chosen 
superposition a |0) + ai|l) of the vacuum and one-photon states. It consists of a three-port 
double-interferometer fed by a one-photon state and a coherent state. The scheme involves 
only linear optical elements and avalanche photodetectors, and therefore should be of interest 
from the point of view of the experimental realization. In principle, i. e. in case of a perfect 
photodetection process, the setup can be used to generate any chosen superposition with 
a conditional probability about 20%. The imperfections of realistic photodetectors have 
been taken into account, and their effects have analyzed in details. An optimal working 
regime has been found, in which output states arbitrarily close to the desired superposition 
is obtained with non vanishing conditional probability. Typical values for the fidelity are 
above F > 99%, with conditional probability about P ~ 1%. For the relevant case of 
balanced superposition, the detection probability may be increased by a fine tuning of the 
amplitude of the coherent input. In this case, the performances of the setup are approaching 
the ideal working regime. 
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APPENDIX A: CONDITIONAL DENSITY MATRIX FOR REALISTIC 

PHOTODETECTION 

Starting from Eqs. (|), (§) and (|20|) we have 
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d n = cos 2 <^|n Y |/?><<5|n N |<5) 

d 00 = sin 2 0i[sin 2 (f) 2 ([3\bfl Y tf\p)(6\fl N \6) + cos 2 2 (/3|rL / |/3)(5|cII N c t |5) 

- sin 2 cos0 2 ((/?|fr Y 6 t |/?)(5|cfl N |(5) + (f3\bfl Y \f3)(6\fl N J\5)^ 
d i = sin 0x cos 0! ( sin 2 {(3 \ fl Y tf \ (3) (5 1 fI N 1 5) + cos 2 (/3 1 ft Y | /3) (5 1 n N c f 1 5) ) ( Al ) 

where (3 = 7cos0 2 and 5 = 7sin0 2 . Let us denote by \z) a generic coherent state with 
complex amplitude z E C, then by using the definition (|20|) of the POM < LTn, fiy r we have 



(z|n N |z) = exp {— r]\z\ 2 } 
(z\afl N \z) = z(l — rj) exp {— ?7|z| 2 } 
;|aII N a^|z) = (1 — rj) l + \z\ 2 {\ — rj) exp {— ^|z| 2 } 



(A2) 



and 



(z\U Y \z) = 1 - (z\U N \z) 
(z\afl Y \z) = z — (z|aII N |z) 
(z|aII Y a^|z) = 1 + |-z| 2 — (z|afl N a^|z) . 



(A3) 



Eventually, upon inserting Eqs. (|A2|) and ( |A3| ) in Eq. 
for the coefficients of the conditional output state. 



we arrive at the expression (|23| 



APPENDIX B: MODELING DETECTORS AS PHOTOCOUNTERS 

By modeling a detector as a photocounter we assume that it is able to discriminate among 
pulses of different amplitudes, ideally corresponding to the different number of recorded 
photons. Actually, the number of " clicks" cannot be not the number of incoming photons, as 
the photocounter is characterized by a nonunit quantum efficiency rj. The POM describing 
the measurement is given by a Bernoulli convolution of the ideal photon-number POM 
fin — \ n )( n \- m formula, we have 
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k—n 



k=n 



\k){k\ 



(Bl) 



Therefore, compared to the picture of detectors as avalanche photodetectors, we have that 
the operator probability for the vacuum detection is the same, i.e. I1 N = ITq, whereas for 
the case of a single click 



Eqs. (|A]|) are now transformed into 





= r]\z\ 2 e~ 


-v\z\ 2 




\z) 


= r]\z\ 2 e~ 


-V\z\ 2 


[l + \ z \ 2 (l-r))] 


\z) 


= r]\z\ 2 e~ 


-V\z\ 2 


[1 + 3|^| 2 (1 — r?) + k| 4 (l 



which leads to a detection probability given by 

PJo = rj sin 2 (pi sin 2 2 + |7| 2 cos 2 2 (1 — rj sin 2 
and to a conditional output state whose coefficients are expressed as 

du = v\l\ 2 cos2 0i cos2 02 

d 00 = T] sin 2 0i sin 2 2 + | T 1 2 ( 1 _ V) cos2 02 

<^oi = ^7 sin 0i sin 02 cos 0i cos 02 1 + 2|7| 2 (1 — rj) cos 2 02 . 

For the balanced setting 0i = 2 = tt/A we have P^ = T]/4[l + |7| 2 (2 — rj)] and 

d n 

doo = "oi 
Finally, the fidelity is given by 



|7| 2 


1 + | 7 | 2 (2 


-v) 


1 + H 2 (1 


-rj) 


1 + | 7 | 2 (2 


-rj) 





1 


4 (1- 


- v) 


(1 + | 7 | 2 )[1 + 


7 


3(2-77)] 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 

which is the result reported in Ref. ||. The fidelity of Eq. ( |B7D varies in the range 5/6 < 
F < 1 as a function of the quantum efficiency at the photodetectors. 
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FIGURES 



a 
b 




FIG. 1. Schematic diagram of the conditional interferometric setup for the preparation of any 
chosen superposition ao|0) + ai|l) of the vacuum and one-photon states. The BS' are identical 
balanced beam splitters, whereas D a and D c denote two identical avalanche photodetectors. The 
first stage consists of a Mach-Zehnder interferometer fed by an one-photon state in mode b. Then, 
one of the output from the first interferometer is used as input of the other one, whose second port 
(mode c) is excited in a weak coherent state. Both the output modes from the second interferometer 
are detected by avalanche photodetectors, and depending on the observed result we obtain different 
conditional output state in the mode a (denoted by OUT in the picture). The event of recording one 
photon in one of the photodetectors (either or D c ) and no photons in the other one corresponds 
to the preparation of a superposition of the vacuum and one-photon states. The relative weights 
of the two components may be tuned by varying either the internal phase-shifts 6\ and 62 or the 
amplitude \ j\ of the coherent input, whereas the the phase of superposition equals the argument 
arg7. 
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FIG. 2. Density plot of the conditional detection probability P* for the preparation of the bal- 
anced superposition of Eq. (|i~5|). The maximum value (P^q ~ 21%) is reached for <fti = (f>2 — 0.715, 
corresponding to an optimum amplitude \"f\ pt — 0.755. However, the dependence of the detection 
probability on (pi and 4>2 is not dramatic, such, as it is apparent from the plot, there exists a 
sizeable region in which the detection probability is above P^ > 20%. 
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FIG. 3. Performances of the setup with realistic description of the photodetectors: the case of 
balanced choice for the internal phase-shifts. The figure shows the detection probability P\q (on 
the left) and fidelity F to the desired superposition (on the right) as a function of the intensity 
of the input coherent state for different values of the quantum efficiency of the photodetectors. In 
both plots we have, from bottom to top, rj = 20%(solid line), 40%, 60%, 80%, and 100%. 




FIG. 4. Performances of the setup with realistic description of the photodetectors: preparation 
of balanced superpositions. The figure shows the detection probability Pio and the fidelity F to 
the desired superposition as a function of the internal phase-shifts <p\ and (j>2 for unit quantum 
efficiency (on the left) and for r\ = 50% (on the right). 
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